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We study the Josephson effect in a superconductor-ferromagnet-superconductor (SFS) junction 
with ferromagnetic domains of noncollinear magnetization. As a model for our study we consider 
a diffusive junction with two ferromagnetic domains along the junction. The superconductor is 
assumed to be close to the critical temperature T c , and the linearized Usadel equations predict a 
sinusoidal current-phase relation. We find analytically the critical current as a function of domain 
lengths and of the angle between the orientations of their magnetizations. As a function of those 
parameters, the junction may undergo transitions between and -k phases. We find that the presence 
of domains reduces the range of junction lengths at which the ir phase is observed. For the junction 
with two domains of the same length, the n phase totally disappears as soon as the misorientation 
angle exceeds £ We further comment on the possible implication of our results for experimentally 
O ' observable 0-7r transitions in SFS junctions. 
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The interest in proximity structures made of superconducting and ferromagnetic layers (respectively, S and F) in 
contact with each other has been recently renewed due their potential applications to spintronics 1 - and to quantum 
computing!^ The interplay between superconductivity (which tends to organize the electron gas in Cooper pairs 
with opposite spins) and ferromagnetism (which tends to align spins and thus to destroy the Cooper pairs) leads to 
a variety of surprising physical effects (for a review, see Ref. 0). As a consequence of the exchange splitting of the 
Fermi level^ the Cooper pair wave function shows damped oscillations in the ferromagnet, leading to the appearance 
of the so-called a ir state" in SFS junctions^ In the ir state, the superconducting order parameter is of opposite sign 
in the two S electrodes of the Josephson junction, and thus a macroscopic superconducting phase difference of tt 
appears in the thermodynamic equilibrium. This phase difference should lead to spontaneous nondissipative currents 
in a Josephson junction with annular geometry.— A possible signature for the 7r-state appearance is a cancellation 
of the Josephson critical current followed by a reversal of its sign as a function of the junction length, - The recent 
experimental observations of critical-current oscillations in experiment o 8 ' 9 ' 10 ^ 11 have demonstrated such 0-ir transitions 
\ as a function of the ferromagnet thickness and temperature. 

The appropriate formalism to deal with mesoscopic S/F junctions has been derived by Eilenberger^ The equations 
of motion for the quasiclassical Green function (averaged over the fast Fermi oscillations) can be further simplified in 
the diffusive regime, i.e., when the motion of the electrons is governed by frequent scattering on impurity atoms: the 
Green functions can then be averaged over the momentum directions. This averaging is justified as long as the elastic 
mean free path l e is much smaller than the relevant length scales of the system, namely the size of the layers, the 
superconducting coherence length £5 = \J D /2nT c , and the length characterizing the Cooper pair wave function decay 
1 . in the ferromagnet £p — y/D/h. Here and in the following, D denotes the diffusion constant, T c the superconducting 
critical temperature, h the magnitude of the exchange field, and the system of units with ^ = fcs=/is = lis chosen. 
The diffusive limit is reached in most of the experimental realizations of S/F heterostructures. In this limit, the Green 
functions can be combined in a 4x4 matrix in the Nambu <S> spin space, and this matrix obeys the Usadel equations? 1 ^ 
SFS Josephson junctions with homogeneous magnetization have been studied in detail within this framework^ Close 
to T c , the proximity effect is weak, the Usadel equations can be linearized, and the current-phase relation is sinusoidal 14 
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Ij = I c sm(j>, (1) 

where <f> is the superconducting phase difference across the junction. The critical current I c shows a damped oscillatory 
dependence on the F-layer thickness (for a review, see Refs. [H and [HI). 

However, understanding the effect of a nonhomogeneous magnetization is of crucial interest for obtaining a good 
quantitative description for the critical-current oscillations in SFS junctions. Indeed it is known that real ferromagnetic 
compounds usually have a complex domain structure. Strong ferromagnets (such as Ni or Fe) consist of domains with 
homogeneous magnetization pointing in different directions whereas the mag netic structure of the weak ferromagnets 
(Cu-Ni and Pd-Ni alloys) used in the experiments reported in Refs. lill^ fTolfllT is still not known precisely. The problem 
of SFS junctions with inhomogeneous magnetization has been addressed previously for spiral magnetizations 1 ^ and 
in the case of domains with antiparallel (AP) magnetizations i 18 d 9 In the latter case, the critical-current oscillations 
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FIG. 1: S-FF'-S junction with non-collinear magnetization. 



(and thus the ir state) are suppressed in the symmetric case where the F layer consists of two domains of the same 
size. This can be explained by a compensation between the phases acquired by the Andreev reflected electrons and 
holes, of opposite spins, in the two domains J£ 

In the present paper, we extend that analysis to the case of a SFF'S junction close to T C} with the two magnetic 
domains F and F' of arbitrary length and relative orientation of the magnetizations. To emphasize the effect of the 
misorientation angle between the magnetizations of the two domains, we choose to minimize the number of parameters 
in the model. The interfaces are then chosen to be perfectly transparent, and spin- flip scattering is neglected in both 
S and F layers. Furthermore, we assume that the diffusive limit is fully reached, that is we do not take into account 
corrections due to a finite mean free path (note that for strong ferromagnets the magnetic coherence length £p may 
become comparable to l e ). 

The main result of our calculation is that, in the symmetric case where the two domains have equal thicknesses, 
we obtain a progressive reduction of the 7r-state region of the phase diagram as the misorientation angle increases. 
Surprisingly, the ir state completely disappears as soon as the misorientation angle 9 exceeds -|. 

The paper is organized as follows. In Section II we solve the linearized Usadel equations and give the general 
expression for the Josephson current. In Section III we discuss the simplest cases of parallel and antiparallel relative 
orientation of magnetizations with different domain sizes d\ and d^. We obtain analytically the full phase diagram 
in d\— d% coordinates. In agreement with Ref. [ill , the 7r state is absent in the symmetric case d\ — di for domains 
with antiparallel magnetizations. In the asymmetric case, the critical current oscillates as di is varied while keeping 
d\ constant. For sufficiently thick layers (<ii,2 3> £f)> the critical-current oscillations behave like in a single domain of 
thickness \di — c?2 1 + (7t/4)£f. In Section IV, we discuss the case of an arbitrary misorientation angle in the symmetric 
configuration d\ = c?2 = d. In the limit when the exchange field is much larger than T c , we derive analytically the 
0-7r phase diagram of the junction depending on the junction length d and and on the misorientation angle 9. We 
show that the ir state disappears completely for 9 > -|. In the last section V, we discuss possible implications of our 
findings for experimentally observed 0-7T transitions in SFS junctions. 



II. MODEL 



We study a diffusive SFF'S Josephson junction with semi-infinite (that is, of thickness much larger than £g) 
superconducting electrodes, as shown in Fig. [TJ The phase difference between the S-layers is denoted <fi — 2\, the 
thicknesses of the two ferromagnetic domains d\ and c?2 • In the following we consider a quasi-one-dimensional geometry 
where the physical quantities do not depend on the in-plane coordinates. For simplicity, we assume that the S-F and 
F-F' interfaces are transparent. We further assume that the temperature is close to T c so that AcT, and this allows 
us to linearize the Usadel equations. 

In the case of superconducting-ferromagnet systems, the proximity effect involves both the singlet and the triplet 
components of the Green's functions^ The Usadel equation in the ferromagnetic layers takes the form (we follow the 
conventions used in Ref. [2lh 

(gVg) - uj [r 3 a , g] - i [f 3 (h • &) ,g] = 0. (2) 

The Green function g is a matrix in the Nambu <g> spin space, f a and a a denote the Pauli matrices respectively in 
Nambu (particle- hole) and spin space, ui = (2n + 1) ttT are the Matsubara frequencies and h is the exchange field in 
the ferromagnet. The Usadel equation is supplemented with the normalization condition for the quasiclassical Green 
function 



g 2 = 1 = tq&o- (3) 
For simplicity, we assume that the superconductors are much less disordered than the ferromagnets, and then we can 
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impose the rigid boundary conditions at the S/F interfaces: 
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Nambu 



where A denotes the superconducting order parameter and the different signs refer respectively to the boundary 
conditions at x = —d\ and x = d 2 . 

Close to the critical temperature T c , the superconducting correlations in the F region are weak^ and we can linearize 
the Usadel equations I©,© around the normal solution g — t^&q. The Green's function then takes the form 
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(5) 



where the scalar f (respectively /J) and vector f (respectively ft) components of the anomalous Green functions 
obey the linear equations 
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The projections of the anomalous Green function on the direction of the exchange field ("parallel" components) are 
defined as f± (x) — /g ± f ^) • where is the unit vector in the direction of the field. The "perpendicular" 

component f_jy refers to the axis orthogonal to the exchange field. Generally, this component is a two-dimensional 
vector. In our system, however, f lies in the plane spanned by the magnetizations in the two domains, and therefore 
/p has only one component. 

It follows from Eq. ^ that the decay of the "parallel" and the "perpendicular" components is governed by two very 
different length scales. The parallel component decays on the length scale while the perpendicular component is 

insensitive to the exchange field and decays on the typically much larger scale £g = C-Fy (experimentally, h may 

be more than 100 times larger than T c , see, e.g., Ref. ITT1 ) . 

In the absence of the exchange field, / CT and /J components are related by complex conjugation. The exchange field 
h breaks this symmetry, and the relation between f a and /J becomes 



The solutions to the equations 



fkx) = U-x)- 

|HJ in each of the ferromagnetic layers are given by 
f± j_(x) = A 3 ± j_ sinh A±,j_x + B± ± coshA±.^x, 



(8) 



(9) 



where the 12 coefficients Ai_ , and B± ± (j 



_i_ aim u ± j_ yj — 1,2 denotes the layer index) must be determined using the boundary 

conditions at each interface. Note that it is enough to solve the equations for the functions p a : the functions fp can 
be then obtained from the symmetry relation ijHJ) 

As we assume transparent S/F interfaces and rigid boundary conditions, 
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= fi( x = d 2 )=0 



(10) 

At the (perfectly transparent) F/F' interface, the standard Kupriyanov-Lukichev boundary conditions 2 - provide the 
continuity relations 



/a(* = 0) = fl{x = $) 



9 J1\ 
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(11) 
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(here a takes values to 3 and refers to a fixed coordinate system). Note that, in the general case, since the 
ferromagnetic exchange fields do not have the same orientation in the two F-layers, the latter conditions do not lead 
to the continuity of the reduced functions /i j_ and their derivative, except in the parallel case. 
The last step will be to compute the Josephson current density using the formula^ 

Ij = ieN(0)DSnT ^Tv (r 3 a gd x g) , (12) 

UJ — — oo 

where S is the cross section of the junction, N(0) is the density of states in the normal metal phase (per one spin 
direction) and the trace has to be taken over Nambu and spin indices. The current can be explicitly rewritten for the 
linearized g 

oo 

Ij = -ieN(0)DSnT ^ 

UJ — — C 

Using the coefficients introduced in equations ©, the Josephson current (fT5|) reads 

Ij = ieN (0)DSttT ]T y [A a { X )B a {- X ) ~ B a { x )A a {- x )\ + X ± [A ± (x)B x (- X ) - B ± { X )A X (- X )} . (14) 

OJ,CT=± 

Since the coefficients A^ and B^ are solutions to the linear system of equations (fit)]) and (fTTj) . they are linear 
combinations of e %x and e~* x . Since the expression (|14[) is explicitly antisymmetric with respect to X i— > — X , it always 
produces the sinusoidal current-phase relation |T]). Finally, the expression (TTJ]) does not contain the domain index j: it 
can be calculated in any of the two domains, and the results must coincide due to the conservation of the supercurrent 
in the Usadel equations. 

In the following sections, this formalism is used to study the influence of a magnetic domain structure on the 
Josephson current. 



\{f*d x fl - fld x f a ) + f ± 8 x f[ - f[d x f ± 



(13) 



III. DOMAINS OF DIFFERENT THICKNESSES IN THE P AND AP CONFIGURATIONS 



A. Parallel case (P) 



In the most trivial case 8 = 0, the equations can be solved easily with A J ± = B 3 ± = 0. We naturally retrieve the 
expression reported in Ref. ,4) for a single-domain S-F-S trilayer (of thickness c?i + ^2), 



ff = eN(0)DSTrT V 

— J / 1) 



U),<7 — ± 



sinh X a {di + d 2 ) 



(15) 



The exact summation over the Matsubara frequencies u> can be done numerically. However, in many experimental 
situations, the exchange field is much larger than T c . In this limit, we can assume /i>w which implies X± = ^Ei. 
The summation over Matsubara frequencies reduces then to 
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and the critical current is given by the simple expression 



7 C P = I Re 
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with 



h = 



(1 + 

eN(0)DSTrA 2 
2^T " 



(17) 



(18) 



From Eq. (jT7J) it is clear that the critical current oscillates as a function of the junction length, with a pseudo-period 
of the order of £p. When the critical current becomes negative, the S-F-S hybrid structure is in the tt state. 
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FIG. 2: Quasi-periodic to ir transitions for antiparallel (solid lines) and parallel (dashed lines) magnetization. On the graph, 
the indications and n refer to the antiparallel case. In the parallel case, the transitions occur along lines with di +^2 = const, 
starting from the state. 



B. Antiparallel case (AP) 

In the antiparallel configuration 9 = n, the exchange field has the opposite direction in the two domains. In this 
case we again find A J ± = = 0, and the critical current can be easily derived, 



2A fT A_ 



A CT sinh \-,?d 2 cosh A CT <ii + \- a cosh \- a d 2 sinh X a d\ 

In the limit of the large exchange field h 3> T c , the summation over the Matsubara frequencies (|16[) results in 

2 



I* p = I Re 



with 



sin(<i + + id_) + sinh(<i + — id-) 



d+ = (di + d 2 )/^ F 
d_ = (di-<h)/fr. 



(19) 



(20) 



(21) 
(22) 



For plotting the 0-7T phase diagram in d\-d 2 coordinates we use the condition of the vanishing critical current. 
From the equations (fT7|) and (|2H)l . the critical current vanishes if 



sin d+ cosh d+ + sinh d+ cos d+ = 



in the parallel case, and if 



sin d + cosh <i_ + sinh d+ cos cL 







(23) 



(24) 



in the antiparallel case. The resulting phase diagram is plotted in Fig. O 

For d\ = e?2 = d (symmetric case), we obtain that the critical current is positive for any d: identical F layers in the 
AP configuration cannot produce the n state (a similar conclusion was drawn in Ref. [TH for ballistic junctions and for 
diffusive junctions at low temperature). For d\ ^ d 2 , the SFF'S junction can be either in the usual state or in the 7r 
state depending on the difference between d-y and d 2 (see Fig. [5]). For large d\ and d 2 , the periodic dependence of the 
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FIG. 3: Critical current dependence on the size of the junction for 8 — 0,7r/4 and 3 7r/4. We take h/T — 100 which corresponds 
to a realistic value for a diluted ferromagnet, as reported in Ref. [ll] and £t = \/D/2irT. 

phase transitions on the layer thicknesses approximately corresponds to a single-layer SFS junction of the thickness 
\di — rfs j + ("7t/4)£f- This result is similar to the case of the clean SFF'S junction where the phase compensation 
arising from the two antiparallel domains is observed^ 

Another interesting feature of the phase diagram in Fig. [2] is the "reentrant" behavior of the phase transition at a 
very small thickness of one of the layers. If the SFS junction is tuned to a 0-7T transition point, and one adds a thin 
layer F' of antiparallel magnetization, then a small region of the "opposite" phase (corresponding to increasing the 
F thickness) appears, before the F-F' compensation mechanism stabilizes the phase corresponding to reducing the F 
thickness. 

In this Section, we have seen that the tt state disappears in the antiparallel orientation for geometrically identical 
F-layers. However, we do not observe an enhancement of the critical current (compared to the zero field current) in the 
AP configuration such as reported in Refs. H9II23I . This in agreement with the claim of Ref. [l^ that this enhancement 
is present only at low temperatures. 

In the next section, we demonstrate that the suppression of the 7r state occurs continuously as we change the 
misorientation angle. 

IV. ARBITRARY MAGNETIZATION MISORIENTATION AND EQUAL THICKNESSES 

In the previous section, we have plotted the phase diagram for arbitrary layer thicknesses d\ and di in the cases of 
parallel and antiparallel magnetization. In principle, one can extend this phase diagram to arbitrary misorientation 
angles 9. Such a calculation amounts to solving a set of linear equations (fTU|) and (fTTj) for the 12 parameters defined 
in Eq. This calculation is straightforward, but cumbersome, and we consider only the simplest situation with 
equal layer thicknesses d\ = di = d. 

For equal layer thicknesses, the 0-ir transitions are present at 9 = and absent at = tt. We will see below 
that with increasing the misorientation angle 9 the amplitude of the critical-current oscillations (as a function of d) 
decreases, and the tt phase progressively shrinks. At a certain "critical" angle 9 C , the tt phase disappears completely 
for any value of d. We find that the critical value is 9 C = | , surprisingly independent of the strength of the exchange 
field. 

The details of the calculation of the critical current are presented in the Appendix. In the general case, the current 
can be written in the form of a Matsubara sum such as given in Eq. (|A5|) . In Fig. [3J we plot the current as a function 
of the domain thickness for different angles performing the summation over Matsubara frequencies numerically using 
realistic values for the temperature and the exchange field. We find that the domain structure reduces the 7r-state 
regions compared to the 9 = parallel case as well as the amplitude of the current in this state. To the contrary, 
the 0-state regions are extended and the current amplitude is increased in this state. This result may be simply 
understood as a continuous interpolation between a sign-changing I c in the single-domain case and an always-positive 
I c in the antiparallel case. 

Considering the high-exchange-field limit introduced in Sec. Ill, namely h T c , and assuming further d <C £5 
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FIG. 4: phase diagram in the limit of a large exchange field. The dependence on d is almost (but not exactly) periodic. 

(which is a reasonable assumption for the first several 0-ir transitions in the high-field limit), we have A^ -C A± and 
X±d <C 1 so that one can expand Eq. IA5|) in powers of Aj_. To the lowest order of expansion, the sum over Matsubara 
frequencies is done and we obtain 



I c (9) = 8^- Jo 



(Q+ + P+ tan 2 §)(P+ + Q+ tan 2 § ) - (1 - tan 4 §)P_ 



£f (P| - P 2 + tan 2 f (P + g + + P_Q_))(Q2 _ Q2_ + tan 2 |( P+ Q + + P_Q_)) ' 

where P± and Q± are simple functions of the ratio rf/^F, 

d / , % c? , . d \ 

P± = 2— cosh(l + i)— ± cosh(l - i) — 

?F V ?F ?F/ 

Q± = (1 + i) sinh(l - i)^- ± (1 - i) sinh(l + i)^-. 

?F ?F 



(25) 



(26) 



From the general formula (f25|) . one can retrieve the expressions (fT7|) and (|20|) for the Joscphson current in the 
(symmetric d± = d%) parallel and antiparallel cases by setting respectively 8 — and 8 = tt. Within the approximation 
of a high exchange field, the critical current (123)) is a ratio of second degree polynomials in the variable tan 2 |. The 
critical current cancels if 



tan 4 -(P+Q+ + P_Q_ 



tan 2 -(P 2 + Ql) + (P+Q+ - P_Q_) 



0. 



(27) 



This equation allows one to compute the full S-FF'-S phase diagram in the d-0 coordinates (Fig. [4]). We observe that 
Eq. (|2T|) cannot be satisfied for any thickness as soon as 8 exceeds As the misorientation angle 9 decreases below 
■|, the region of the tt state in the phase diagram Fig. [4] increases, and it becomes maximal at 8 — (i.e., in the 
parallel configuration). 

Away from the high-exchange- field limit, we can find the value of 8 C numerically using the exact formula, Eq. (|A5|) . 
Remarkably, our calculations show that the critical value 8 C = f remains independent of the strength of the field h. 



V. DISCUSSION AND EXPERIMENTAL ASPECTS 



The main conclusion of the present work is that a domain structure in the SFS junction reduces the region in the 
phase space occupied by the tt state. We have demonstrated this reduction with the example of the two domains 
placed along the junction. However we expect that this qualitative conclusion survives for more general configurations 
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of domains. In view of this reasoning, we suggest that the shift in the 0-7r transition sequence reported in Ref. [ll| and 
attributed to a "dead layer" of the fcrromagnet may be alternatively explained by a domain structure in the junction. 
If such a domain structure slightly reduces the region of the it phase in favor of the phase, this would shift the 
positions of the two first 0-7T-0 transitions in a manner similar to the effect of a "dead layer" (see, e.g., our 8 = ir/4 
plot in Fig. [3]). To distinguish between the two scenarios, one would need to observe at least three consecutive 0-7r 
transitions (and/or develop a more realistic theory of the effect of domains in SFS junctions). 

Many of our results are based on the high-exchange-field approximation assuming £s/£f = \fh] (27rT c ) 3> 1. This is 
a reasonable approximation for the type of samples reported in Ref. I 111 : the exchange field in the CuNi ferromagnetic 
alloys has been estimated at 850K, whereas the critical temperature of Nb is of the order of 9K. Thus, the ratio £s/€f 
is of the order of 4. Note that the high-field limit is consistent with the diffusive limit condition hr e <C 1 with r e the 
elastic mean free time (see Ref. [2J for details). The parameters of the experiments^ yield the estimation hr e « 0.1. 

In our treatment we have neglected the finite transparency of the interfaces, the finite mean free path of electrons, 
the spin-flip and spin-orbit scattering. Of course, those effects may be incorporated in the formalism of Usadel 
equations in the usual way (see, e.g., Refs. I25ll26ll27t ). We expect that they do not change the qualitative conclusion 
about the reduction of the tt state by the domain structure. A realistic quantitative theory of SFS junctions may 
need to take those effects into account, in addition to a more realistic domain configuration in the ferromagnet. 
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APPENDIX A: SOLVING THE LINEARIZED USADEL EQUATIONS 



To solve the system of linear equations (fTT)|) and (fTTj) with the 12 variables A J ± ± and B± ± , it is convenient first to 
reduce the number of variables by resolving the continuity relations (jilt in terms of the 6 variables 

8 8 
f3± = B\_ =F B]_ tan - = B\±B\ tan - 

, B\-B\ 8 , Bl-B 2 8 
p± = B\ + tan- = B\ - tan- 

8 8 
a± = X±A± qp X±A]_ tan - = \±A± ± X^A\ tan - (Al) 

, Al X+Al-X^Al 8 x a2 X+A 2 -X-A 2 _ 8 
a± = X±A ± + ± tan - = X±A]_ ^ tan -. 

Solving now the set of 6 equations (fTT))) produces the solution 



a±_ = — A + A_Aj_cosh(A ± d) — — =- tan -(1 + tan 2 -)cosx 



\ \ \ . .1. , \ , ... , 1 _|_ ^ ar * 

/3+- 0- = -— A + A_sinh(A_L<2) — -(1 + tan 2 -) cosx 



V\ -P- +tan 2 %(p+q+ +p~q~) 2 2 



4A p_ 

— A + A_sinh(A_L<2) — —„ > 

w p + — P- + tan z ~(p + q + + p-q-) <t 



4A P+ + tan 2 f<7+ 

/3++/3_ = — A+A_ sinh(A ± ci) — ^— ^ rcosx 

uj p z + — pt + tan z %(p+q+ + P-q~) 



4iA 



A + A_ Aj_ cosh (Xj_d) — ^- -(1 + tan 2 -)sinx (A2) 



q\ - q- +tan 2 ^{p+q+ + P-q-) ' ' " 2' 
2 e„ 



4iA q++tan z |p + 

= A + A_A ± cosh (Aj_a) — 5-3- smx 

w q\ - q z _ + tan 2 %{P+q+ + P-q-) 



where 



2iA. A , /A q- 6 , 2 0. 

A+A_ smb. (X±d) — — — — 5-3- ■ tan -(1 + tan -) sin X ■ 

ui q+-q 2 -+ tan 2 |(p + q+ + P-q-) 2 2 



p± = A + A_ sinhAj_(i(coshA+(i± coshA_d) 

q± = X± cosh X±d(X + sinh X-d± X- sinh X + d) . (A3) 
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In terms of the new variables, the supercurrent (|14p becomes 

1 



I j = ieN(0)DSnT ^ 

UJ 

1 



-(a++a-)(x)09++/3-)(-x) 



+ 



4(1 + tan 2 f ) 



(o + -a_)(x)09 + -/9_)(-x) + 



(Xx(x)P±(-x) 



1 + tan 2 | 



be -x] ■ 



The resulting current-phase relation is sinusoidal with the critical current 

(A + A_) 2 A^ sinh2A ± d 



I c (0) = Uot F T 2 J2 



(A4) 



(A5) 



(q+ +p + tan 2 §)(f>+ + 9+ tan 2 | ) - (1 - tan 4 § )p_<?_ 
(p 2 , -p 2 , +tan 2 f - <? 2 + tan 2 §(p+g + 

Eq. (|A5[) can be used for numerical calculations of the critical current for an arbitrary relative orientation of the 
ferromagnetic exchange fields and for any value of their magnitude (e.g., for producing the plot in Fig. [3]). In the 
body of the article, a simpler expression for the current is given in the high-exchange-field limit (Eq. (|25[) ). 
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